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Abstract: We use the AdS/CFT correspondence to study AA = 4 supersymmetric SU{Nc) 
Yang-Mills theory, in the limits of large Nc and large 't Hooft coupling, coupled to a number 
Nf of massless hypermultiplet fields in the fundamental representation of the gauge group. We 
identify a [/(I) subgroup of the R-symmetry under which the fermions in the hypermultiplet 
are charged but the scalars are not. All the hypermultiplet fields are also charged under 
a U{1) baryon number symmetry. We introduce an external magnetic field for the baryon 
number U{1), which triggers the spontaneous breaking of the U{1) R-symmetry, and we 
then introduce a chemical potential for the U{1) R-charge, producing a state with a nonzero 
density of the U{1) R-charge. The system should then exhibit superconductivity of the U{1) 
R-charge. The dual supergravity description is a number Nf of D7-branes in AdS^ x 
with angular momentum on the and a worldvolume magnetic field. We study the zero- 
temperature thermodynamics of the system, and find that for sufficiently large magnetic 
field the system prefers to be in the symmetry-broken phase. For smaller magnetic fields we 
find a discontinuous free energy, indicating that our gravitational setup does not capture all 
equilibrium states of the field theory. 
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1. Introduction 

The anti- de Sitter / Conformal Field Theory (AdS/CFT) correspondence [1-3], and more 
generally gauge-gravity duality, provides a novel tool for studying strongly-coupled systems. 
In particular, gauge-gravity duality provides tractable examples of strongly-coupled gauge 
theories, and thus may shed light on the low-energy dynamics of the theory of the strong in- 
teractions. Quantum Chromodynamics (QCD). More recently, gauge-gravity duality has been 
used to study physical phenomena familiar from solid state physics, such as the quantum Hall 
effect [4,5] and superconductivity [6-12]. Gauge-gravity duality may thus provide solvable 
toy models for condensed matter systems. 

In this paper we will: 1.) identify a specific gravitational system whose field theory dual 
includes a global U{1) symmetry under which only fermions are charged, 2.) study states 
of that theory with a finite density of the U{1) charge, and 3.) study a mechanism that 
triggers the spontaneous breaking of the U{1) symmetry, such that the system should be in 
a superconducting state. We will not study the transport properties of the system, however: 
in this paper we will study only zero-temperature, finite-density thermodynamics. 

We will study a conformal field theory, namely = 4 supersymmetric SU{Nc) Yang-Mills 
theory (SYM), in the 't Hooft limit ^ oo with the 't Hooft coupling A = gyjyfNc fixed, 
and in the additional limit A 1. We will introduce a number Nj of AA = 2 supersymmetric 
hypermultiplets that transform in the fundamental representation of the gauge group, i.e. 
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flavor fields. An M = 2 hypermultiplet contains two complex scalars and two Weyl fermions 
of opposite chirality. We will call the scalars "squarks" and the fermions "quarks," in analogy 
with (supersymmetric) QCD. We will work in the "probe limit," in which we keep Nf fixed 
as we send Nc ^ oo, so that Nf Nc, and work to leading order in the small parameter 

Nf/N,. 

The AdS/CFT correspondence is the statement that J\f = 4: SYM theory in the lim- 
its above is equivalent to type IIB supergravity on the near-horizon geometry of very many 
D3-branes, which is AdS^ x S^, where AdS^ is (4-|-l)-dimensional anti- de Sitter space and 
is a five-sphere [1-3]. The correspondence is "holographic" in the sense that the (3-1-1)- 
dimensional field theory is dual to gravity in the (4-|-l)-dimensional AdS space. The Nj 
hypermultiplets appear in the supergravity description as a number Nf of D7-branes embed- 
ded in the AdS^ x geometry [13]. 

A/" = 4 SYM theory has an SO{6)ji R-symmetry (hence the subscript). The flavor fields 
explicitly break this to SO{4)ji x U{l)ii. Of the fields in the hypermultiplet, only the fermions 
are charged under the U{1)r. Indeed, the U{1)b. acts as a chiral symmetry, rotating left- and 
right-handed quarks oppositely, just like the U{1) axial symmetry of QCD. We may thus 
introduce a chemical potential for the f/(l)i? and produce a state with a finite density of 
quarks only. 

To break the U{1)r spontaneously, we will exploit the fact that the hypermultiplet fields 
have a second U{1) symmetry. With Nf flavors of mass-degenerate hypermultiplets, the SYM 
theory will have a global U{Nf) symmetry. We may identify the overall, diagonal U{1) as 
baryon number, U{1)b- We may introduce a non-dynamical, background magnetic field, B, 
for the U {1)b- Previous AdS/CFT studies have shown that such a background magnetic field 
causes spontaneous breaking of the U{1)r symmetry [14-17]. We thus want to study a state 
with a finite U{l)ji density of quarks and a U{1)b magnetic field. 

What is the gravitational description of such a state? The SYM R-symmetry is dual 
to isometries of the S^, and a finite R-charge density in the SYM theory is dual to a state 
in supergravity with angular momentum in the directions [18-21]. To introduce a U{1)r 
density of quarks, then, we must study a D7-branc rotating in the directions. D7- and 
D5-branes spinning in AdS^ backgrounds have been studied in refs. [22,23], respectively. The 
U{1)b symmetry is dual to the U{1) gauge invariance on the worldvolume of the D7-branes, 
hence to introduce a magnetic field in the SYM theory we must introduce a magnetic field 
on the worldvolume of the D7-branes. 

Notice that the U{1)r symmetry will not be broken due to finite density physics, rather, 
we will break the symmetry by an external mechanism, the U{1)b magnetic field. Contrast 
this with the usual physics of Bose-Einstein condensation as described by a theory of a single 
complex scalar field with a potential including a mass term and a quartic term. The theory 
has a global U{1) symmetry which shifts the phase of the field. Roughly speaking, a chemical 
potential for the U{1) acts as a negative mass-squared. If the mass-squared is positive, so 
that the potential has a single minimum and classically the field has zero expectation value, 
a sufficiently large chemical potential causes a second-order phase transition to a state of 
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broken symmetry: the potential changes to a "wine bottle" shape. Our system is analogous 
to a scalar field with a negative mass-squared, in the sense that the symmetry is broken 
even at zero chemical potential. In particular, if our system does describe a superconductor, 
any "pairing" mechanism will be due to the magnetic field, and not the due to high-density 
physics. 

In addition, as the U{l)ji is an R-symmetry, some of the fields of the = 4 multiplet 
also carry the charge, and hence may also contribute to a state of finite U{1)r density. 

We will ignore such a contribution. In supergravity language, we will only study D7-branes 
spinning in AdS^ x S^, rather than D7-branes spinning in the near-horizon geometry of 
spinning D3-branes [18-21]. In SYM theory language, we will introduce different values of 
the chemical potential for different fields: zero for the adjoint fields and nonzero for the flavor 
flelds. We discuss the relevant background (of spinning D3-branes and = 4 SYM theory 
with U{l)ii chemical potentials [24-26]) below, in section |2.3| . In particular, the Af = 4 SYM 
theory (in flat space) with a U{1)r chemical potential has no genuine equilibrium ground 
state [24-26], so by ignoring the chemical potential for the adjoint flelds we are ignoring a 
known instability of the theory. 

Despite its shortcomings, we hope that our system may serve as a nice toy model, i.e. as a 
laboratory for questions about gravitational descriptions of superconductivity. This system is 
attractive mainly because it is relatively simple and because the dual field theory is explicitly 
known (in particular, we know that in the flavor sector only fermions carry the U (1)_r charge). 
Furthermore, we believe this system is more attractive than some other systems with similar 
features, for a number of reasons. 

Other gravitational systems dual to fleld theories with only fermions in the flavor sector 
include for example the D4/D6 and D4/D8 models of refs. [27,28]. An important difference 
between our system and these systems is the potential between heavy test charges in the 
field theory. The D4/D6 and D4/D8 systems describe field theories with confining potentials, 
which are of course preferable when the goal is to study QCD. The AA = 4 SYM theory with 
massless N = 2 flavor fields in the probe limit is conformal, hence the potential between heavy 
test charges is necessarily Coulombic, which may be preferable for some condensed matter 
applications. More generally, condensed matter systems near quantum critical points are 
typically described by strongly-coupled conformal field theories (see ref. [29] and references 
therein). 

Other gravitational descriptions of field theories with spontaneously broken U{1) chiral 
symmetry include for example the D4/D6 system as well as D7-branes in the Constable- Myers 
background [30]. The benefit of using the U{\)b magnetic field to break the symmetry is 
that, in some sense, it is intrinsic to the D7-brane, i.e. in field theory language we can change 
the scale of chiral symmetry breaking (which is determined by the magnetic field B) without 
changing any other physics. In contrast, consider the D4/D6 system of ref. [27]. In that case 
the scale of chiral symmetry breaking and the scale of confinement are both set by the same 
parameter, the Kaluza-Klein compactification scale. If we want to change the scale of chiral 
symmetry breaking, then, we must also change the glueball masses. 
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Our main results are that for sufficiently large U{1)b magnetic field the system prefers to 
be in a symmetry-broken phase. For smaller values of the magnetic field we find a gap in the 
free energy, indicating that our supergravity setup is missing something. More specifically, 
our ansatz for the D7-brane embedding does not seem to capture all values of the free energy. 

This paper is organized as follows. In section ^, we describe the SYM theory and its 
supergravity dual in more detail. In section |3| we study D7-branes rotating in directions 
with no U{1)b magnetic field, extending the analysis of ref. [23]. In section ^ we review the 
results of refs. [14-17] for non-rotating D7-branes in AdS/CFT with worldvolume magnetic 
fields. In section |5| we study the full problem of rotating D7-branes with worldvolume magnetic 
fields. We conclude in section ^ with some discussion and with suggestions for future research. 
We collect some technical details in the Appendix. 

2. The Theory and Its Supergravity Dual 
2.1 The Theory 

We will study the maximally supersymmetric Yang-Mills theory in (3-1-1) dimensions, = 4 
SYM theory, with gauge group SU{Nc). The fields of the AA = 4 supermultiplet include the 
gluons, four Weyl fermions and three complex scalars. The AA = 4 SYM theory is conformal, 
so that the 't Hooft coupling, A = gy^^Ac, is a free parameter. We will take the 't Hooft limit 
of Ac ^ oo with A fixed, followed by the strong-coupling limit A ^ 1. 

We will also introduce a number Nj of AA = 2 supersymmetric hypermultiplets in the 
fundamental representation of the gauge group, i.e. flavor flelds. In the language of AA=1 
supersymmetry, the N =2 hypermultiplet contains two chiral multiplets of opposite chirality: 
two Weyl fermions of opposite chirality and two complex scalars. We will colloquially refer 
to the fermions as "quarks," and the scalars as "squarks," in loose analogy with QCD. 

We will keep Nf fixed as we take Nc ^ oo, and work to leading order in the small pa- 
rameter Nf/Nc- This is known as the probe limit. In the language of perturbation theory, we 
are discarding diagrams that contain quark or squark loops. More physically, we are ignoring 
quantum effects due to the flavor flelds because such effects are parametrically suppressed by 
powers of Nf/Nc- In particular, we ignore the running of the coupling: the beta function of 
the theory is proportional to ^"^jf, which vanishes in our limit, for fixed A. 

The AA = 4 SYM theory has an SO{6)r R-symmetry. The fiavor fields explicitly break 
this to SO{4:)r X SO{2)r. We will denote the current associated with the S0{2)r = U{1)r 
as J'^. What are the charges of the fields under the C/(1)r? A table of charge assignments 
appears in many places, for example in refs. [31,32]. In the AA = 4 multiplet, one complex 
adjoint scalar has charge +2. Two of the adjoint fermions have charge -|-1, and two have 
charge —1. In the fiavor sector, the squarks are neutral under the [/(l)ij. The fermion in 
one N =1 chiral multiplet has charge -|-1, while the fermion in the other chiral multiplet has 
charge —1. 

The U{\)r thus mimics the U{1) axial symmetry of QCD, and we will refer to it as a 
chiral symmetry. Notice in particular that a quark mass term will explicitly break the f/(l)_R 
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symmetry. Additionally, if Nf is on the order of then the U{1)r is anomalous, just like 
the ^7(1) axial symmetry of QCD. In the probe limit, however, the anomaly is not apparent, 
for the same reasons that the running of the gauge coupling is not apparent: the quantum 
effects that produce the anomaly are paramctrically suppressed by powers of Nf/Nc- 

If the quarks are massless, so that the U{1)r is a symmetry of the Lagrangian and J'^ is 
conserved, d/^J^^ = 0, then we may introduce a chemical potential for the U{l)ji. In a state 
with a finite U{l)ii density, any contribution that the flavor fields make to the density must 
come from the quarks, and not from the squarks. Notice that, as the U{1)r symmetry is 
an axial symmetry, a state with a net ^7(1)^ density of quarks is a state with an excess of 
left-handed quarks (for example). 

2.2 The Supergravity Dual 

We begin with type IIB string theory, where we consider an intersection of Nc coincident 
D3-branes and Nf coincident D7-branes described by the array 

Xq Xi X2 Xs X4 X5 Xq X7 Xg Xg 

D3 X X X X (2.1) 
D7XXXXXXXX 

We first consider the D3-branes alone. We take the near-horizon limit of the D3-brane ge- 
ometry, which gives us AdS5 x S^, with radius of curvature L given by = AirggNca'^, 
where gg is the string coupling and a' is the square of the string length. We take the usual 
limits of A^c ^ 00 with ggNc fixed, followed by ggN^ S> 1, so that ^ a'^. In particular, 
in the latter limit, heavy string modes decouple, and we may approximate the string theory 
by type IIB supergravity in AdSr^ x S^. The AdS/CFT correspondence is then the statement 
that supergravity on this background is equivalent to the low-energy theory on the D3-brane 
worldvolume, which is A/" = 4 SYM theory in the 't Hooft limit with large 't Hooft coupling. 

The Nf D7-branes introduce additional open string degrees of freedom, producing fields 
in the fundamental representation of the SU (A^c) gauge group on the D3-branes' worldvolume. 
If we keep Nf finite as we take Nc — >■ oo, so that Nf <^ N^ we may neglect the D7-branes' 
contribution to the stress-energy tensor in the supergravity theory^. The D7-branes thus do 
not deform the geometry: they are probes embedded in AdSs x S^. The D7-branes will be 
extended in the AdS^ directions as well as along an C S^. 

We will use an AdS^ x metric suited to the symmetries of the D7-branes, 

ds"^ = -| rj^^dx^dx" + — (dr^ + r'^ds% + dy^ + y'^dcf)'^) (2.2) 

L Tg 

Here re is the distance to the D3-branes in the transverse M^, dels is the metric of a unit- 
radius S^, and dy^ + y^d(f)^ is the metric of the xg-xg plane written in polar coordinates. 

'^Additionally, the D7-branes source the dilaton and axion, but again, in the probe limit we neglect this 

effect. This is dual to the SYM theory statement that when Nf -C Nc the quantum effects of the flavor fields 
that cause the running of the coupling and the U{1)r anomaly are suppressed. 
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Notice that rg = + y^. The boundary of AdS^ is located at rg — > oo. Starting now, we wiU 
use units in which L = 1 unless stated otherwise. We then translate between string theory 
and SYM theory quantities using a'~^ = AirggNc = Qym^c = A. 

The part of the D7-brane action that will be relevant here is the Dirac-Born-Infeld (DBI) 
term, 

Sd7 = -NfTor J d\^-det {gj^J + i2TTa')Fab) (2.3) 

Here Tdy is the D7-brane tension, are the worldvolume coordinates, gj^J is the induced 
worldvolume metric, and Fab is the U{1) worldvolume field strength. 

The D7-brane has two worlvolume scalars, y and (f). An ansatz for the D7-brane scalars 
that preserves the Lorentz invariance of the Minkowski space directions, and the 5*0(4) x 
50(2) isometry, is Fab = 0, (p = and y{r). The induced D7-brane metric is then 

dslj = rl r^^^dx^'dx" + ^ [dr^ (1 + y'{rf) + ds%) , (2.4) 



and the D7-brane action is 



Sd7 = -MV^z,i j drr^'^/l + ^/W (2.5) 



where we have defined the constant 

(2vr) 



M = NfTmVsz = jT^NfN, (2.6) 



where V53 = 27r^ is the volume of a unit-radius 5^ and in the second equality we have 
converted to field theory quantities using T07 = (2iif ~ Starting now, we will drop 



the factor V^3.i from eq. ( |2.5| ), which represents the (infinite) volume of (3+l)-dimensional 
Minkowski space, and re-define Sm as an action density. 
The equation of motion for y(r) is 

which restricts the asymptotic form of solutions to be 

2/(0 = ^0 + 5 + O (1) (2.8) 

with constant coefficients cq and C2- Here cq is the asymptotic separation between the D3- 
branes and the D7-branes in the -xg plane. A string stretched between the D3-branes and 
D7-branes, whose endpoint represents an excitation in the fundamental representation on the 
D3-brane worldvolume, will have minimum length cq. We may thus identify the mass m of 
the hypermultiplet fields as this length times the string tension: m = = ^ cq. 

The field y{r) is dual to an operator Om in the SYM theory given by taking of the 
SYM theory Lagrangian. The operator Om thus includes the mass operator of the quarks, m 
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times the mass operator of the squarks, and a cubic couphng between the squarks and one 
complex scalar of the = 4 multiplet (the scalar with charge +2 under the U{l)ji). Om 
is written explicitly in the Appendix. We also show in the Appendix that the sub-leading 



coefficient in eq. (p^, C2, is related to the expectation value of Om 



{Om) = -j^^NfN,2c2 (2.9) 

Notice that Om is charged under the U{l)fi symmetry (just as the quark mass operator 
is charged under the U{1) axial symmetry of QCD), and hence, when cq = 0, a nonzero 
expectation value for Om signals the spontaneous breaking of the U{1)r symmetry. 

The D7-brane action, eq. ( |2.5[) , depends only on the derivative of y{r), and hence the 
quantity in parentheses in eq. ( |2.7D is a constant of motion. Solutions with nonzero values of 
the constant of motion have been studied in refs. [33,34]. These solutions are not supersym- 
metric, and in fact describe D7/anti-D7 configurations. For more details, see refs. [33,34]. 

The factor under the square root in eq. ( |2.5| ) is a sum of squares, so the solution with 
the smallest value of the on-shell action, which must be the physically preferred solution, is 
y'{r) = 0, or y{r) = cq, in which case the constant of motion is zero. These embeddings are 
supersymmetric. In the SYM theory, they describe hypermultiplet fields with mass m = ^ cq 
and a vacuum state in which {Om) = 0. 

In the supergravity picture, the y{r) = cq solutions describe D7-branes that "end" some- 
where in AdS^, which is most easily seen from the induced D7-brane metric evaluated on such 
a solution, 

dslj = (r2 + eg) r]^^dx>'dx'' + , (dr^ + ds%) ■ (2.10) 

(r"^ + Cq) 

At the boundary, r oo, the induced metric approaches AdS^ x S^. When r ^ 0, however, 
we see that rg = r'^ -|- ^ Cq , so when cq is nonzero, the D7-brane does not reach rg = 0, 
and we say that the D7-brane ends at rg = cq. From the induced metric eq. ( |2.1C1| ), we can 
see that when r = 0, the S*^ has zero volume. What is happening as r decreases is that the 
C "slips" or "shrinks," as allowed by topology, and eventually collapses to a point at 
r = 0. Notice that if cq = (dual to massless hypermultiplets), then the does not slip, 
and the D7-brane is present at all values of rg. 

The Ricci scalar R{r) associated with the D7-brane's induced metric in eq. ( 2.1C| ) is 



At the boundary r — > og, R{r) —14, which is indeed the curvature of AdS^ x 5'^. At 
the endpoint r ^ 0, R{r) —8. In later sections we will compute the Ricci scalar of the 
D7-brane numerically, and these limits will provide useful checks. 

In subsequent sections we will encounter embeddings for which the curvature diverges 
at r = 0. For a general solution y(r), an easy way to understand such a singularity is by 
looking at the induced metric in eq. (p.4|). To avoid an angular deficit, and hence a conical 



-7- 



singularity, at r = 0, we must have (1 + y'{r)'^) 1, hence y'{r) 0. Solutions for which 
y'(0) is nonzero are thus singular, as we will see below. When the curvature grows we expect 
curvature corrections to the DBI action eq. to become important, therefore we must 

discard such embeddings as unphysical: they are solutions to an equation of motion that arises 
from an action that is no longer a reliable approximation to the actual D7-brane action. 

We want to study a state of the SYM theory with a nonzero U{1)r density. The global 
U{1)r symmetry of the SYM theory is dual to the SO{2) isometry that rotates xg and xg 
into one another, or equivalently that shifts 4>hy a constant. Notice for instance that a finite 
Co explicitly breaks the 5*0 (2) isometry in the x^-xg plane, which is dual to the statement in 
the SYM theory that a finite mass m explicitly breaks the U{1)r symmetry. A state of the 
SYM theory with finite U{l)ji charge density is dual to a supergravity solution with nonzero 
angular momentum in <j). We thus want to study a D7-brane spinning in the xg -xg plane, 
with a non-trivial y{r): if y{r) = 0, so that the D7-brane sits at the origin of the xg-xg plane 
for all r, then obviously the D7-brane will have zero angular momentum in (/>. Moreover, in 
the SYM theory we want massless hypermultiplet fields, so that the U{1)r is a symmetry of 
the Lagrangian, hence we want solutions with cq = 0. For such solutions, the leading term in 
the asymptotic form of y{r) will be the C2/r'^ term. 

To summarize: our goal is to find embeddings in which the D7-brane rotates in 4> and has 
zero Co with nonzero C2, i.e. the D7-brane has zero asymptotic separation from the D3-branes 
but nonzero angular momentum in the -xg plane. 

Notice what supergravity is telling us about the SYM theory: for a D7-brane with cq = 
to have nonzero angular momentum, C2 must be nonzero. Translating to SYM theory 
language, in any state with a finite U {1)r density, the U{l)ji symmetry must be spontaneously 
broken, as indicated by a nonzero (Om). 

In AdS/CFT, we identify the supergravity action, when evaluated on a solution, with the 
generating functional of the SYM theory [1-3]. More specifically, with a nonzero chemical 
potential we identify the on-shell supergravity action Sd7 with the thermodynamic potential 
in the grand canonical ensemble, as Sdj = For spinning D-branes, we identify the 
angular frequency of rotation, uj, with the chemical potential, fi, and the angular momentum 
with the density, (J*) [18-21]. In SYM theory language, the latter is given by (J*) = — 
so in supergravity language, we have (J*) = '^^^^ . Notice that we will thus be studying 
densities proportional to the factor M in Sd7, that is, we will be studying densities of order 
(J*) oc XNfNc- We write an explicit formula for in the Appendix. 

Embeddings for probe D5-branes spinning in AdS^ x were studied in ref. [23]. In 
that case, the dual SYM theory includes flavor fields confined to a (2+l)-dimensional defect. 
The Lagrangian of this theory is written explicitly in refs. [35,36]. Again, in that theory, 
of the fields in the fundamental representation, only the fermions carry the relevant U{l)ji 
charge. Our results for the D7-brane in AdS^ x in section ^ will be similar to those for 
the D5-brane. Our results in section || for the D7-brane with a worldvolume magnetic field 
will be new. 
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2.3 The Adjoint Contribution to the Density 



As mentioned in section 2A_, one of the complex scalars and all of the fermions of the = 4 
multiplet are charged under the U{1)r, hence these fields may contribute to a state with a 
finite U{l)ii density. A correct thermodynamic analysis must include all microstates that 
produce the same macroscopic charge density, hence a correct thermodynamic analysis must 
include states in which the fields of the = 4 multiplet contribute to the density. In other 
words, in the grand canonical ensemble we choose values of the temperature and chemical 
potential, and the dynamics of the theory then determines the ground state. 

In supergravity language, including the adjoint fields means allowing the D3-branes to 
spin in the xs -xg plane. The near-horizon geometry of spinning D3-branes is known [18- 
21]. A complete supergravity analysis would thus involve allowing D7-branes to spin in the 
background produced by spinning D3-branes. Notice, however, that we may give the D3- 
branes and D7-branes distinct angular frequencies. If the frequencies are equal, then in the 
SYM theory we have a single chemical potential for the C/(l)ij. If the frequencies are distinct, 
we have different values of the chemical potential for different fields in the theory, one value for 
the fields of the = 4 multiplet and another value for the fields of the N = 2 hypermutliplet. 
Translating the SYM theory thermodynamic analysis to supergravity language: a correct 
analysis would mean studying a system of D3-branes and D7-branes spinning with the same 
angular frequency and finding the solution that extremizes the on-shell supergravity action. 

As done in ref. [23], however, we will simply ignore the rotation of the D3-branes. In 
SYM theory language, we will introduce a nonzero f/(l)_R chemical potential for the flavor 
fields only. From a SYM theory point of view, then, what we will do is artificial: we introduce 
a f/(l)_R chemical potential only for the quarks, and then use the U{\)b magnetic field to 
force them to pair and thus break the C/(l)ij. Nevertheless, as mentioned in the introduction, 
we hope that this system may serve as a toy model for answering questions about holographic 
superconductors. 

To place our analysis in context, and to understand what ignoring the U{l)ii chemical 
potential for the adjoint fields really means, we will now briefiy review the results of refs. 
[24-26], where the = 4 SYM theory in the presence of a C/(l)/j chemical potential was 
studied at both weak and strong coupling^. One result of refs. [24-26] was that, in fact, for 
any finite U{1)r chemical potential the SYM theory has no equilibrium ground state. 

First consider the AA = 4 SYM theory in the large-A^c limit, in fiat space, at zero tem- 
perature, and at zero 't Hooft coupling. A U{l)ii chemical potential will act as a negative 
mass-squared for the scalar charged under the U{1)b,- In the presence of a C/(1)r chemical po- 
tential, then, the potential has no minimum; the theory has no equilibrium ground state. At 
finite coupling, the superpotential has a moduli space parameterized by mutually-commuting 



^Notice that refs. [24-26] focus primarily on the TV = 4 SYM theory formlated on a spatial three-sphere, so 
that the adjoint scalars acquire a curvature coupling that acts as a positive mass-squared. The phase structure 
then becomes more interesting than for the theory in flat space: for the details, see the phase diagrams in the 
references. Roughly speaking, we can obtain the phase diagram for the theory in flat space by taking a "large 
volume" limit in which the radius of the three-sphere goes to infinity (relative to all other scales). 



- 9 - 



constant background values for the adjoint scalars, and indeed, the zero-temperature behav- 
ior persists from weak to strong couphng [25,26], i.e. the theory at zero temperature has no 
equiUbrium ground state. 

At finite temperature, a weak-coupUng analysis of the A/" = 4 SYM theory with a U{1)r 
chemical potential has been performed in refs. [24,26]. The principal result was that for 
chemical potentials below a critical value /icrit = 

VXT, the origin of the moduli space is 
meta-stable. More precisely, when < ficrit, the potential exhibits runaway behavior for 
large values of the scalars, but the origin of the moduli space is a local minimum with a 
lifetime that grows exponentially with Nc- The meta-stability was discovered by computing a 
one-loop effective potential for the scalars, plotted against the expectation values of the scalar 
eigenvalues. The potential barrier between the meta-stable state and the unstable state is 
lowest in the case of a single eigenvalue splitting from the rest. For /j, > iJ,critj the potential 
barrier disappears, and with it the meta-stable state. 

The finite-temperature story is qualitatively the same at strong coupling, where the 
system can be analyzed using AdS/CFT. In particular, in ref. [25], an analysis of probe D3- 
branes spinning in the near-horizon geometry of spinning D3-brancs revealed that the meta- 
stability persists to strong coupling: roughly speaking, the stack of Nc spinning D3-branes 
"spits out" individual D3-branes one at a time. The supergravity picture thus nicely agrees 
with the field theory picture of a single eigenvalue separating from the rest and penetrating 
the potential barrier. 

To return to our system: introducing D7-branes in the probe limit will not alter the 
physics of the background produced by the spinning D3-branes since in the probe limit we 

ignore the back-reaction of the D7-branes. If we included the rotation of the D3-branes, then, 
we know a priori that the system is cither meta-stable or unstable: it cannot be the ground 
state because the system has no ground state. In other words, if we did include the rotation of 
the D3-branes, then we would expect the D3-branes to carry most of the angular momentum, 
and indeed to exhibit runaway behavior^. 

D7-brane probes in the near-horizon geometry of spinning D3-branes have been studied 
in refs. [22,37]. The principal result was that the nonzero U{1)r chemical potential does not 
trigger spontaneous breaking of the U{1)b. symmetry^. The case of D7-branes with worldvol- 
ume magnetic fields probing the near-horizon geometry of spinning D3-branes remains to be 
studied. 

With the above background in mind (and in particular, remembering what we are ignor- 
ing), we now turn to our analysis of D7-branes spinning in AdS^ x S^. 

^We will mention in passing that we can "fix" the instability: in SYM theory language, we can compactify 
the spatial directions into a three-sphere, which, for the J\f = 4 SYM theory alone, stabilizes the theory for 
sufficiently small chemical potential. 

''Notice that the result of ref. [37] seems to be the opposite: that the U{1)r chemical potential causes 
breaking of the U{1)r in the flavor sector. As indicated in refs. [22,23], however, that conclusion came from 
using unphysical D7-brane embeddings. In fact, a U{1)r chemical potential does not cause breaking of the 
U{1)r. 
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3. Finite ^7(1)/? Chemical Potential 



To study D7-branes spinning in AdS^ x S^, we consider the following ansatz for the D7-brane 
worldvolume scalars: y(r) and (t){t,r) = uot + f{r). We will thus have a D7-brane spinning 
with frequency oj in the xs-xg plane. As explained in ref. [23], and similar to the system 
in ref. [38], the r dependence in <j3{t,r) is required to guarantee the reality of the D7-brane 
action for all values of r, for certain embeddings. 

We can also motivate the r dependence in <l){t,r) via T-duality [22]. If we perform a 
T-duality in the (p direction, the D7-brane becomes a D8-brane and (f){t,r) — > Affj{t,r), hence 
the D8-brane now has a constant electric field pointing in the <j) direction: Ff 



t<f> 



UJ. 



From 



previous experience with electric fields on D-branes in AdS/CFT [38], we expect that, to 
guarantee reality of the D7-brane action for all values of r, A^i, must have radial dependence 
of the form A(jj{t,r) = ut -\- f[r) and hence, T-dualizing back to the D7-brane, we find the 
(t){t, r) written above. 

With our ansatz, the induced metric g^^ of the D7-brane has components 



D7 11 
grr = 9rr + OyyV + 



n 



9tt + 



12 



yrt 



(3.1) 



where primes denote differentiation with respect to r and dots denote differentiation with 
respect to t, and with all other components identical to eq. (|2.4|) . The D7-brane action 
becomes 



Sd7 = -N / (irr%/(l + y'2) ( i_<^2__J^_^ ) +y2^/2 



(y2 +r2)2 



We will define the Lagrangian C via Sdi = — f dr C (notice the sign). The action depends 
only on (p'{r), so the system has a constant of motion, which we call c. 



_6£_ 



(l + 2/'2) (^l_a;2^-^j+y20'2 



c. 



(3.3) 



We then solve algebraically for (j)'{r), 



(p'{r) = - 

y 



J\J'2y2^6 — 



(3.4) 



Plugging this into the action, we find 



Sd7 



-M / dr y^TT y 



12 

\ 



1 _ 

1 _ c2 i 



(3.5) 
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We can derive the equation of motion for y{r) either by varying the action in eq. (^]^) and 
then inserting the solution for (p'ir), or by varying the Legendre-transformed action Sd^, 



Sd7 = Sd7 - I dr(j)'(r) ^^^'^ 



J dr 



64>'{r) 



y2 \ / 



We will not write the equation of motion explicitly. 

The numerator and denominator under the square root in the action eq. ( |3.5| ) can change 
sign as r goes from infinity to zero. If one changes sign while the other does not, then the 
action will become imaginary. Both must change sign simultaneously for the action to remain 
real, i.e. the numerator and denominator under the square root in eq. ( |3.5D must share a 
common zero. We thus find two curves in the (r, y) plane, and the D7-brane must either cross 
both simultaneously or cross neither for the action to remain real. The first curve is (here we 
restore factors of the AdS radius L) 

1-^ (^2+>)2 =0 (3-7) 
Which is just the equation for a semicircle of radius ^wL^, centered at (0, ^ujL"^): 

y-^u^L^y+r' = \iv'L\ (3.8) 

As noted in ref. [23], the entire D7-brane spins with constant angular velocity uj, and its linear 
velocity ujy{r) depends on r, as does the local speed of light, which decreases as r decreases. 
The semicircle is where the D7-brane's linear velocity equals the local speed of light. The 
second curve is a cubic, 

2/(r) = ^^. (3.9) 

Notice that the semicircle is determined only by the value of so once we choose to it is the 
same for all solutions. The value of c, however, varies from solution to solution. For example, 
we will show shortly that some solutions reach r = without ever crossing the semicircle, so 
the numerator under the square root in eq. ( |3.5| ) remains positive for all r. For the action to 
remain real, the denominator under the square root must also be positive for all r, which is 
only possible if c = 0. For solutions that do cross the semicircle, the value of c is fixed entirely 
by the position where the D7-brane crosses, that is, if (ro, yo) is the point where the D7-brane 
crosses the semicircle (so that tq and yo obey eq. (^)), then c = N yor^. In other words, 
every solution that crosses the semicircle has its own cubic curve. The qualitative behavior of 
c as a function of the y position on the semicircle is depicted in fig. |l|. The maximum value 
of c occurs at y = I ujL? , and c goes to zero at y = and y = ujL"^ . 
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From the equation of motion, we find the 



asymptotic form of y{r), 



y{r) 



Co + 



£2 



+ 77 Co + O 

2 



log r 



(3.10) 

Notice that when cq is nonzero, a finite lo pro- 
duces a logarithmic term at order l/r"^. As 
mentioned above, we are interested in solu- 
tions with Co = 0, for which the logarithmic 
term will be absent. We show in the Appendix 
that, given a solution for y(r) with nonzero w, 
the expectation value of Om is determined by 
y(r)'s asymptotic coefficients as 




Figure 1: The value oic/N as a function of the 
y-position on the semicircle of eq. (3.S). Here we 
have chosen a; = 1 as an example. 



(a 



(^2c2 + ^w'co + ^ ^0 log(c^)^ . (3.11) 



From the explicit solution for (j)'{r), we can also find </)(t,r)'s asymptotic form, 



ujt 



1 c 1 
2 



+ O 



log r 



(3.12) 



Notice the factor of Cq in the denominator of the coefficient of the 1/r^ term, which suggests 
that solutions with co = (the ones we want) must have c = 0, as otherwise (j){t, r) diverges 
asymptotically. Indeed, we have found numerically that this is always the case. 

The field (j) is dual to a SYM theory operator that we will denote O,^, which is the phase 
of the hypermultiplet mass operator Om- We write explicitly in the Appendix. We also 
show in the Appendix that the constant c determines the expectation value of Off, as {Off,) = c. 

Our numerical solutions for y{r) are depicted in fig. |2|. We generate these as follows. We 
divide solutions into two classes, those that reach r = "above the semicircle," and those 
that intersect the semicircle. 

Solutions that end above the semicircle have c = 0. These solutions describe D7-branes 
for which the collapses to zero volume before intersecting the semicircle. To generate 
these, we specify the value of y{r) at r = 0, which must be > uj, and we require that y'{r) 
vanish at r = 0, to avoid a conical singularity, as explained in section 2.2. 



Solutions that intersect the semicircle have nonzero c. The equation of motion depends 
on c, so to generate solutions we first need to fix c, which we do simply by choosing a point 
(ro,yo) on the semicircle. We then need the value of the derivative at the semicircle, y'(ro). 
As in ref. [23], we can derive a regularity condition on y'{r) from the equation of motion 
expanded about the semicircle. We have found numerically, however, that the solutions are 
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amost entirely insensitive^ to the value of y'iro). Given (ro, yo) and y'{ro), we can numerically 
solve the equation of motion for all r. 

The solutions in fig. ^ are qualitatively 
similar to the solutions found for spinning 
D5-branes in ref. [23]. The solutions are 
also qualitatively similar to the solutions 
found in refs. [17,39] for D7-branes with a 
worldvolume electric field, dual to the SYM 
theory in the presence of external U{1)b 
electric field, which is easy to understand 
via T-duality arguments such as the one we 
gave above. We will make three comments 
about the solutions. 

First, as found in ref. [23] for spinning 
D5-branes, the only solution with cq = is 
the trivial solution, y{r) = 0, which has no 
angular momentum. In other words, using 
SYM theory language, for massless quarks, 
introducing the U{l)ji chemical potential does not break the U{1)r, but also does not produce 
a U{l)ji density. We have not found a good SYM theory argument for why this is so. Notice 
also that all the solutions that cross the semicircle have nonzero cq. 

Second, all the solutions with nonzero cq describe flavor fields in the SYM theory with 
time-dependent masses. More specifically, as (p corresponds to Otp, the phase of the mass 
operator Om, these solutions describe time-dependent masses of the form me*"^*. We do not 
have a good field theory intuition for the physics of such a mass term. 

Third, all of the solutions that cross the semicircle are singular at r = and hence should 
be discarded as unphysical. We can see the singularity easily from fig. |2|: these solutions all 
have nonzero y'{r) at r = 0. Additionally, given our numerical solutions we have computed 
the Ricci scalar associated with the induced metric, eq. and observed the divergence 

explicitly. 

In fig. ^ we illustrate the behavior of the Ricci scalar for solutions ending above the 
semicircle. As r — > oo, we see R{r) — > —14 for all solutions, the expected value for AdS^ x S^. 
Solutions that end far above the semicircle, with y(0) ^ uj, should approach the constant 
solution y(r) = cq of a non-spinning D7-brane, and hence at r = should have R{0) — > —8. In 
fig. ^ we see that is the case. As y{0) decreases toward lo, however, we see that the curvature 
at r = decreases, and appears to diverge when the D7-brane ends precisely at the semicircle, 
y(0) = Lo. Such behavior is in fact familiar [30,40-46]: probe D7-branes in AdS-Schwarzschild 
may end "above" the black hole horizon or may intersect the horizon. The "critical solution" 
that ends precisely at the horizon is singular. We are seeing the same behavior, with the 
black hole horizon replaced by the semicircle. 

^In more detail: despite changing the value of y'{r) at the semicircle by five orders of magnitude, and even 




Figure 2: Numerically-generated solutions for 
spinning D7-brane embedding functions y{r). The 



semicircle of eq. (3.8) is also depicted. Here we have 
chosen w = 1 as an example. 



- 14 - 




(a) (b) 

Figure 3: The behavior of the Ricci scalar R{r) for D7-branes that end above the semicircle of eq. 
(^), i.e. for which j/(0) > oj. We have used w = 1 to generate these figures, (a.) The Ricci scalar for 
solutions with various values of 2/(0), ranging from y(0) = + 15 (the top curve) to y(0) — w + 0.76 
(the bottom curve). All solutions have R{r) —14 as r — > oo. (b.) The Ricci scalar evaluated at 
r = 0, i?(0), as a function of the position above the semicircle where the D7-brane ends, y(0). We 
see that when y{0) ^ uj, R{0) w —8 as appropriate for the constant solution y{r) ~ cq. As y{0) 
approaches the top of the semicircle, y{Q) uj, we see that i?(0) diverges to negative infinity. 



We will not present plots for the Ricci scalar of solutions that intersect the semicircle. We 
will only note that, again, as r ^ oo, all solutions have R{r) —14, and that the curvature of 
every solution diverges at r = 0, as expected. Indeed, the curvature appears to be extremely 
large everywhere inside the semicircle. For example, setting co = 1, the solution intersecting 
the semicircle at the point (rg, yo) ~ (0.29, 0.91) reaches a curvature on the order of 10^ inside 
the semicircle within a distance 0.2 of the semicircle, and a curvature of order 10'' at a distance 
of 0.4. Clearly curvature corrections to the DBI action will be important for such solutions, 
so we cannot trust them. Nevertheless, we will include such solutions in our later analysis, 
for two reasons. First, these solutions are required to account for the full range of cq values. 
Second, we expect that, in an AdS-Schwarzschild background, the high-curvature region may 
be cloaked by the horizon, in which case such solutions may become physically acceptable. 
The analogous figures in refs. [17,23,39] for D-branes in AdS-Schwarzschild suggest this. 

As discussed in refs. [22,23], this system undergoes a first-order phase transition in which, 
roughly speaking, as y{0) approaches to the D7-brane "jumps" from ending outside the semi- 
circle to intersecting the semicircle. The transition is analogous to the D7-brane's first-order 
"meson melting" phase transition in the AdS-Schwarzschild background [30, 40-47] (for de- 
tails, see refs. [22, 23]). The transition is between two solutions with nonzero cq. We are 
interested only in solutions with cq = 0, so we will not investigate the phase transition here. 



changing its sign, the solutions always "settle down" to the solutions depicted in fig. g within a very short 
distance from the semicircle. 
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4. Finite U{1)b Magnetic Field 



With Nf massless flavor fields, the SYM theory has a global U{Nf) symmetry. We identify the 
overall diagonal U{1) as baryon number, U{\)b- In the supergravity description, the U{1)b 
current is represented by the U (1) gauge field, A^, propagating on the D7-brane worldvolume. 
We can describe external electric and magnetic fields in the field theory, coupled to anything 
carrying U{\)b charge, by introducing non-normalizable modes for in the supergravity 
theory. For example, we will be interested in a magnetic field, which we introduce by adding 
to our D7-brane ansatz the constant field strength F^y = B. In the SYM theory, we identify 
Fxy as a constant U{1)b magnetic field pointing in the z direction. 

The utility of introducing B is that, at zero temperature, zero mass, and zero U{\)b 
chemical potential, AdS/CFT calculations have shown that the U{\)b magnetic field triggers 
spontaneous breaking of the U{l)ii symmetry [14-17]. In supergravity language, the D7-brane 
is "repelled" from the origin of the (r, y) plane, so that the solution with zero asymptotic 
separation, cq = 0, is no longer just the trivial solution y{r) = 0. 

To illustrate how this occurs, we will briefly review the results of refs. [14-17]. We consider 
an ansatz for the D7-brane fields with y{r) and Fxy = B only (so for now (f){t,r) = 0). We 
will also define the notation B = (2Tra')B. The D7-brane action is then 



= _^ dr 



„3 



From the equation of motion, we find the asymptotic form of y{r), 

y{r) = CO + 5 + O (1) , (4.2) 

where again we translate to SYM theory quantities with m = and (Cm) oc — 2c2. 

We generate solutions numerically as follows. For all solutions we impose y'(0) = 0. We 
then choose the value of y(0) and numerically integrate to large r. From these solutions we 
extract the values of cq and C2. 

We present the plot of C2 as a function of cq in fig. |^ (a.). The curve actually spirals into 
the origin, crossing the vertical axis an infinite number of times. We thus have infinitely many 
solutions with cq = 0. As argued in ref. [17], however, the cq = solution with lowest energy 
will be the physical one, which turns out to be the "first" cq = solution, "first" meaning the 
first Co = solution we reach as we enter the plot from the right (from large values of cq). In 
fact, the other co = solutions are not only thermodynamically disfavored, they are unstable, 
having tachyonic fluctuations [15]. Notice in particular that the trivial solution y(r) = 0, at 
the center of the spiral, is unstable. 

The physical cq = solution has nonzero C2, indicating that in the SYM theory the 
[/(l)ij is spontaneously broken. An analysis of the D7-brane's fluctuation spectrum, dual to 
the SYM theory's meson spectrum, confirmed the existence of a Goldstone boson associated 
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yfr) 

0.6 r 




(a) (b) 

Figure 4: (a.) The value of the sub-leading asymptotic coefficient C2 as a function of the leading 
coefficient cq for a D7-brane with nonzero worldvolume magnetic field B = [2'Ka')B. We have used 
B = 1. The curve actually spirals into the origin, and hence crosses the cq — Q (vertical) axis an 
infinite number of times, although that is not apparent in our plot. The physical cq = solution has 
C2 = 0.226. b.) The solution y{r) for the physical cq = solution, for _B = 1. 



with the symmetry breaking [14]. In the SYM theory, when m = the only scale in the 
problem is B, hence by dimensional analysis we have {Om) oc B^/"^ . More precisely, (Cm) = 
- (2^ ^iV/A^c (2 X 0.226) ^3/2 [14]. 

A picture of the physical cq = solution in the (r, y) plane appears in fig. ^ (b.). This is 
just what we want: a solution with zero asymptotic separation but nonzero extension into the 
(r, y) plane. Our goal, roughly speaking, is to set this solution spinning, giving us a solution 
with zero asymptotic separation but nonzero angular momentum. 

5. Finite U{1)r Chemical Potential and U{1)b Magnetic Field 

We will now study D7-branes spinning with angular frequency w, and with a constant world- 
volume magnetic field F^y = B, which will produce embeddings with cq = but C2 7^ 0. 
Such solutions describe states in the SYM theory with massless hypermultiplet fields, a finite 
U{1)r density, and spontaneous breaking of the U{1)r symmetry. 

We consider an ansatz for the D7-brane worldvolume fields with y{r), (j){t, r) = ujt + f{r) 
and now F^y = B. The DBI action becomes 



The constant of motion c is now 

.r^^ yV(l+(pf^, ^^^^^ 
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The solution for (f)'(r) is now 



y 



(5.3) 



— c 



Plugging this into the action gives 



Sd7 



AA / dr rM 1 + 



^y2 _|_ 

The Legendre transform of the action is 



,12 



\ 



,j2 y 



1 + 



52 



(5.4) 



M'^ y'^ 



s 



D7 



S 



D7 



dr (j)'{r) 



SSd7 

6(1)' (r) 



-M 



.12 

\ 




2^6 



(5.5) 



We can see that the equation for the semicircle is unchanged, but the cubic curve has 
become 

/ ^2 1 d2 

(5.6) 




y(0 = Vtto ifi 



The value of c is still fixed uniquely by a point on the semicircle. The plot of c versus y is 
qualitatively similar to fig. |l|. 

The asymptotic forms of y{r) and (j){t,r) are unchanged from the B = case, eqs. ( ^.10 ) 
and ( 3.12| ). (Cm) and (O^) are again given by eq. ( p. 11 ) and {O^p) = c, respectively. 

We generate solutions numerically in precisely the same way as in section ^. We first 
consider solutions that intersect the semicircle. The behavior with nonzero B is more com- 
plicated than with zero B, so in fig. ^ we present only a few examples. We choose two points 
on the semicircle and generate solutions with increasing B. As B increases, for the solution 
intersecting the semicircle near the top, the value of cq first decreases, but then begins to in- 
crease. For the solution intersecting the semicircle near the bottom, the value of cq decreases, 
but only very little. We summarize the behavior of solutions with the three-dimensional plot 
in fig. ^, where we plot cq as a function of B and the value of y where the solution intersects 
the semicircle, which we denote yo- 

The behavior of the Ricci scalar for the semicircle-intersecting solutions is qualitatively 
similar to the B = case: for all the solutions that intersect the semicircle, R{r) diverges at 
r = 0, and indeed is extremely large almost everywhere inside the semicircle, and hence the 
solutions must be discarded as unphysical. 

The principal result of figs. ^ and ^ is that all solutions that intersect the semicircle, 
which have nonzero c, also have nonzero cq. We are only interested in solutions with cq = 0, 
however, so we will pay no more attention to solutions that intersect the semicircle. 
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Figure 5: Illustration of behavior of semicircle-intersecting solutions as B increases. Here we set 
w = 1. The top figure has i? 0, and the subsequent figures show the behavior as B increases. We 
present plots for the values B — 3,6.73 (the critical value), 9 and 24. 

Now for the solutions that end above the semicircle. We present some examples of such 
solutions in fig. ^, for increasing values of B. Here we find that for a given point above the 
semicircle the corresponding value of cq decreases as B increases. Heuristically, as we increase 
-B, the solutions "bend down." Turning things around, if we imagine fixing cq and integrating 
into the bulk, then as we increase B we see that the point (above the semicircle) where the 
solution reaches r = increases (y(0) increases). 

At a critical value of i? ?s 6.73, a cq = solution appears: see fig. |^. A cq = solution 
continues to exist for larger values of i?, as we show in fig. ^ (a.). Indeed, if we look only 
at the Co = solution, and increase B, we find that the solution's value of y(0) increases. 
Heuristically, the D7-brane "bends out" more into the y direction as B grows. Notice that, 
from the SYM theory point of view, when cq = the only scales in the problem are B and 
so when working with cq = solutions we will always write B in units of u). The critical 
value of B is thus B w 6.73 w^. 

The trivial solution y{r) = is another cq = solution, so once B juP' reaches the critical 
value we can make a meaningful comparison between two cq = solutions. To determine 
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Figure 6: Three-dimensional plot of co as a function of B and the value of y where the solution 
intersects the semicircle, denoted uq. Here we use uj ~ 1. The only points of the surface that reach 
Co — are at yo — 0, corresponding to the trivial solution (which intersects the semicircle at yo = 0). 

which is preferred, we must, in supergravity language, compare the values of their on-shell 
action Sdt, or, in SYM theory language, compare the values of their free energy 0. Recall 
that we identify Sd7 = so that the solution with larger So7 will be thermodynamically 
preferred. 

The on-shell action Soi suffers from divergences coming from the integration over the 
infinite volume of AdS^. In SYM theory language, these are UV divergences, which we can 
cancel with counterterms. In the Appendix we perform the "holographic renormalization" of 
the on-shell action by regulating and then cancelling the divergences using counterterms. We 
will denote the renormalized action Sren-, and identify Sren = 

The trivial solution has Sren = 0. We find numerically that when B/uJ^ = 6.73, the 
non-trivial cq = solution has Sren/-^ ~ 11-3 > 0, and hence the nontrivial solution is 
thermodynamically preferred. The value of Sren for the nontrivial cq = solution increases 
monotonically with B/uP' as shown in fig. ^ (b.), so the non-trivial solution remains the 
preferred solution as we increase B/uP' . 

Naively, then, we might think the system exhibits a first-order phase transition. In 
supergravity language, the D7-brane "jumps" from the trivial embedding to a nontrivial 
embedding, with nonzero angular momentum. In SYM theory language, the theory jumps 
from a state in which {0„l) = {J^) = to a state with nonzero (Om) and (J*). In fig. ^ (c.) 
we plot (Om) (divided by (27ra')AA) as a function of B/lo'^ and in fig. ^ (d.) we plot (J*) 
divided by as a function of B /uP , both for solutions with cq = 0. Clearly both are nonzero 
at the critical value B juP' = 6.73. As {Om) and (J*) are first derivatives of we seem to 
have a first-order transition. 
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Figure 7: Illustration of the behavior of solutions that reach r = above the semicircle. Here we 
have set lli ~ 1. For a given value of j/(0), as we increase B the associated value of cq decreases. We 
present plots for the values B = 0,3, 6.73 and 9. At the critical magnetic field B « 6.73 the first 
solution with cq = appears. 

That is not correct, however, because the free energy itself is discontinuous: it jumps 
from zero to nonzero at the critical value of BjiiP' . Such behavior is unphysical, and signals 
to us that something is missing. More specifically, some class of cq = embeddings appears 
to be absent for values oiBjiJ^ below the critical value. What kind of embeddings could "fill 
the gap" is not obvious to us, so we leave this as an open question. 

Moreover, we should not trust the nontrivial embedding precisely at the critical value of 
B/uj^ because such an embedding has very high curvature. Indeed, the story of the scalar 
curvature for embeddings that end above the circle is qualitatively the same as in the B = Q 
case: the curvature is finite everywhere, but the curvature at r = grows as y(0) approaches 
uj (the top of the circle), where it diverges. In fig. ^ (a.) we plot the Ricci scalar for the 
Co = solutions as we increase the value oi B/uP' . For the critical value B juP' , the curvature 
of r = diverges to negative infinity. In fig. |9| (b.) we plot the value of the Ricci scalar at 
the r = endpoint for the cq = solutions as a function of B/u}"^. We see again that the 
curvature diverges at r = for the critical solution and then increases as we increase B/uP. 
The closer we come to the critical solution, the less we should trust our solutions. 
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y(r) 




(c.) (d.) 

Figure 8: (a.) The behavior of solutions with cq = as BjuJ^ increases. We use values of BjuP' 
ranging from B/uP' = 6.73 (bottom curve) to B /uP = 30 (top curve). As B /uJ^ grows, the associated 
value of ?/(0) increases: the D7-brane "bends out" more in the y direction, (b.) The value of the 
renormalized D7-brane action divided by A/" as a function of BjuP' for solutions with cq = 0. The 
value of Sren/J^ at B/uj"^ — 6.73 is approximately 11.3. (c.) The expectation value {Om), divided by 
{2TTa')Af, as a function of BjuP' for solutions with cq = 0. (d.) The density (J*) divided by A/" as a 
function oi BjJ^ for solutions with cq = 0. 



6. Conclusion 



We have numerically constructed solutions, reliable within the supergravity approximation, 
for a spinning D7-brane with a worldvolume magnetic field embedded in AdS^ x . These 
solutions describe = 4 SYM theory coupled to massless M = 1 hypermultiplets in a 
state with a nonzero background \J(X)b magnetic field and nonzero ?7(1)_r charge. The 
C/(l)s magnetic field causes spontaneous breaking of the \J(\)r symmetry, hence the system 
should exhibit \J(X)r superconductivity. We initiated the study of the zero-temperature 
thermodynamics of the system, and determined that for large enough values of the magnetic 
field the system prefers a state of broken symmetry. For smaller values we found that our 
class of D7-brane embeddings was insufficient to describe all equilibrium states of the SYM 
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Figure 9: The behavior of the Ricci scalar for D7-brane solutions with finite uj and B and with 
Co = 0. (a.) The Ricci scalar R{r) for the co = solutions for values of B/uj^ ranging from BjiiP' = 6.73 
(the bottom curve) up to B joJ^ = 30 (the top curve), (b.) The value of the Ricci scalar at the r = 
endpoint, R{0), for Co = solutions as a function of B/ui'^. 

theory. We will end with some suggestions for future research^. 

The biggest open question is, of course, where are the "missing" cq = embeddings 
for small values of B/uj'^, which must fill the "gap" in the free energy that we discovered in 
section I? We have left this as an open problem, being content that we found solutions with 
the properties we wanted (cq = 0, nonzero angular momentum, and finite curvature). 

An important generalization would be to introduce a finite temperature T, corresponding 
to D7-branes spinning in an AdS-Schwarzschild black hole background [48]. We have hope 
that many of the singular embeddings we found would be "cured," in the sense that the 
high-curvature region would fall behind the black hole horizon. Additionally, in refs. [16,17], 
the finite-temperature physics of D7-branes with zero lo but nonzero B was studied, with 
the result in the SYM theory that at high temperature the U{l)ji symmetry is restored. We 
expect something similar to happen when we include nonzero oj. In SYM theory language, 
such a transition should extinguish the U{l)ji superconductivity. 

Similarly, another important generalization would be to introduce a background magnetic 
field for the U{1)fi, symmetry. We expect on general grounds that a sufficiently large magnetic 
field should extinguish superconductivity, that is, should restore the U{1)r symmetry. More 
generally, the phase diagram in the plane of U{1)r magnetic field versus temperature should 
be explored. 

The fluctuation spectrum of spinning D7-branes should be computed, which should re- 
veal the effects of a finite U{l)ii chemical potential on the meson spectrum of the SYM 

®In our suggestions we continue to ignore the chemical potential of the adjoint fields. Of course a worthwhile 
extension would be to include the adjoint fields' chemical potential, and to study the system on a spatial 
three-sphere, which stabilizes the theory for sufficiently small chemical potential [24-26]. A good question, 
in supergravity language, is whether our spinning D7-brane with worldvolume magnetic field minimizes the 
on-shell action, as opposed to a state in which the background geometry carries the angular momentum. 
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theory. Such an analysis should also exhibit explicitly the Goldstone boson associated with 
the breaking of the U{l)ii symmetry. 

We chose an ansatz for the worldvolume fields that preserved many symmetries, such 
as translation invariance in the field theory directions. An interesting extension would be 
to consider a more general ansatz, respecting fewer symmetries. Indeed, in QCD at low 
temperature, asymptotically high baryon number chemical potential, and in the large-A^c 
limit (with Nf fixed), the ground state may break translation invariance, forming a so-called 
"chiral density wave" [49,50]. 

We have claimed that our system describes a superconductor, so perhaps the most exciting 
task for the future would be studying the transport properties associated with the U{l)ii 
charge and exhibiting superconductivity explicitly. We expect, for example, to see a gap in 
the frequency dependence of the U (1)r conductivity. 
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Appendix: Holographic Dictionary 

In this appendix we will: 1.) write explicitly the SYM theory operators dual to the D7-branc 
worldvolume fields y and (p, 2.) regulate and renormalize the on-shell D7-brane action, and 
3.) compute the expectation values of the SYM theory operators dual to y and (p, as well as 
the expectation value (J*), from the renormalized D7-brane action. 

First, we will identify the operators dual to y and cp. These have been identified in several 
references; we borrow the results of ref. [51]. We decompose the M = 2 hypermultiplet into 
two A/" = 1 chiral multiplets of opposite chirality. Let Tp and q denote the Weyl fermion and 
complex scalar of one chiral multiplet, and -0 and q the Weyl fermion and complex scalar of 
the other chiral mutliplet. In particular, ^ and ip have opposite chirality. 

The operator dual to y is the supersymmetric completion of the mass operator, which we 
denote Om- In terms of the SYM theory fields, is 

Ora = ii>i^ + q (m + V2^^ q^ + q (m + V2^^ q^ + h.c. (1) 

Here we use the notation $ to denote the complex scalar of the J\f = 4 multiplet with U{1)r 
charge +2. 

The field (p is dual to the phase of Om {i-e. fluctuations of (p arc dual to fluctuations 
of the phase of the hypermultiplet mass term). We denote this operator as 0<^. In terms of 
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SYM theory fields, is 



11)11) + iV2q^ + iV2q^ + h.c. 



(2) 



Next, we will show how to compute finite on-shell D7-brane actions. The on-shell action 
diverges due to integration over the radial coordinate, r all the way to the AdS^ boundary 
at r = cxD. To obtain a finite on-shell action, we first regulate the integral by introducing a 
cutoff, r = A. We then introduce counterterms localized at the r = A hypersurface to cancel 
the divergences of the action, and then send A — > oo, yielding a finite result. This procedure 
is called "holographic renormalization" [52-57]. We will denote the induced metric on the 
r = A hypersurface as j^^, 



dsl^f^ = 7^1, dx^dx'^ = A^ rj^^dx^dx'' 



(3) 



and it determinant as simply 7, so that ^—7 = A^. 
We will give (j) the general coordinate dependence 



(x, r) = • x + /(r) 



(4) 



with x^ the coordinates of (3+l)-dimensional Minkowski space (hence n runs from to 3) and 
a four-vector: = (—a;, h) with spatial vector k. We will denote purely spatial indices 
with lower-case Latin indices, for example, a component of k will be /cj. We will also use the 
notation ((9(/>)^ = r]'^'^ d^(p du(p- 

We denote the regulated D7-brane action as Sreg, so that Sreg = — dr C. With the 
ansatz y{r), (j){x,r) and F^y = B, we have 



-'reg 



M j drr^ 



\ 



+ y^'^ 1 + 



_|_ j>2\2 



(5) 



Inserting the asymptotic forms of the solutions in eqs. (|3l0|) and (with uj"^ ^ - {d(pf), 

we find 



^reg 



A 

M I dr 
1 



1 1 ^ol 



+ -4 id<pr- + -B'- + o 



-N 



-k^ + -cl {d(t)f logA^ + i^^logA^ + O 



r 2 r 

1 
4 



log r 



log A 



(6) 



We then introduce counterterms on the r = A hypersurface to cancel the divergences. The 
counterterms that we use throughout this paper are 



L2 = -^^fV^ (7'^"9^</'9.</.) 500(A) (log 500(A) + 1) 



(7a) 
(7b) 
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1 



^3 = +-N ^/^{2TTa'Y F'^Fij log(A 



(7c) 



where g^(f){K) = and Fij is (27ra') times the D7-brane worldvolume field strength, which 
for us will include only Fxy = B. Written more explicitly, the counterterms are 



L2 = -\^rcl {d^f (log(cg) - log(A^) +1) + o (i^) 



+-^MB'' log(A2) 



where the divergences in A cancel those in the regulated action, eq. 
on-shell action, Sren-, is 



). The renormalized 



(9) 



A number of finite counterterms are possible. Indeed, some of the terms in eqs. (0) 
are finite. We have chosen the particular counterterms above so that {Om) will have the 
appropriate behavior as cq — > oo: in SYM theory language, when m — > oo the flavor fields 
decouple from the dynamics of the SYM theory. We must have {Om) ^ in this limit, which 
fixes the finite counterterms. We will ignore all other possible finite counterterms (i.e. we 
will set their coefficients to zero). 

We will now compute the expectation values (Om) and (O^j,). In the AdS/CFT corre- 
spondence [1-^3], we identify the on-shell supergravity action with the generating functional 
(or grand canonical potential) of the SYM theory as Sren = We thus have 



{O, 



5n 

6m 



-{2'KOi 



I 6Sre 



'6y{A) 



The contribution from Sreg is 



6S, 



reg 



Sy{A) 



r=A 



2c2 + - {dcpf CO (1 - log(A2)) + O 



log A 



(10) 



(11) 



where in the second equality we have inserted the asymptotic solutions eqs. (|3.10| ) and ( [3.121) . 
Of the counterterms, Li and L3 contribute nothing, while the contribution from L2 is 

'2/(A)' 



6L2 
6yiA) 



1 



N{d(t>Yy{A) 



log 



A2 



+ 2 



\n {d^f CO [log(c2) - log(A2) + 2] + O 
Summing everything and taking A ^ 00, we find 

{(Dm) = -(27ra')AA (2C2 - ^ id<pf cq - ^ (dcjyf cq log(c2) 



(12) 
(13) 

(14) 
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In terms of SYM theory quantities, the prefactor is {2'Ka')M = ^fXNf N^. 

To show that (Om) ^ as cq — > oo, we need to know the large-co behavior of C2. To 
determine this, we borrow arguments of ref. [14]: to study large cq, we let y{r) = cq + Y{r) 
and linearize the Y{r) equation of motion, retaining only the leading terms in (cg + r^)"^. 
The equation of motion for the fluctuation Y{r) is then 



2co^2 + CO {dcpf {cl 



which has the solution 



y(r) = ai 



02 



4 



(eg + r2)3 

2 log(cg + r2 



Cn + 







1 B^o 



4 r2(cQ + r^) 



(15) 



(16) 



with integration constants ai and 02 • To fix ai we demand that lim^^oo y{r) = (so that 
Y(r) does not alter the value of cq), which fixes ai = 0. To fix 02, we argue that for sufficiently 
large cq, the solution for Y{r) must be valid for all r. In particular, Y{r) must be finite as 
r — > 0, which means that the coefficient of the l/r^ term must vanish as r ^ 0, which means 
we must have 

1 .9 1 .9 . , 9, 1 ^2 



as = - CO (90)^ + ^co (90)2 log(c2) + - — 
2 4 4 Co 

Inserting this value of 02 into y(r) and then extracting C2 from the r 

C2 = ^ CO {d^f + \co {d<t>f \0g{cl) + \^. 

00 we indeed have (Cm) - 



(17) 

00 limit, we find 

(18) 

0, for our particular 



We can then see from eq. (|14D that as co 
choice of counterterms. 

For the cq = solutions that we want the finite counterterms, and the cq terms in (Cm), 
will not contribute. In numerical calculations, however, we necessarily deal with solutions for 
which Cq is nonzero. We have used the counterterms above for all of our numerical calculations. 
Now for (O^), which is simpler. We take 

¥(A) 



¥(A) 



lim 

A— >oo 



E 



5Li 



\ _ 5£ 

5(t){K)) ~ W 



c. 



(19) 



r=A 



where in the final equality we have used eq. (|3.3D for the constant of motion, = c. Notice 
that the counterterms contribute nothing. 

Lastly, we want to compute the C/(1)_r density (J*). In the SYM theory we have (J*) = 



so in supergravity language we have (J*) 



dSr 



du) 



We will now restore (90)^ 



-ijj 



Let us first compute the contribution from S'reg, borrowing arguments from ref. [16]. We start 
with the regulated action, eq. (|5|), evaluated on a solution. The action is a functional of the 
fields y{r) and 4)(t,r), and when evaluated on a solution has explicit to dependence as well as 
implicit dependence through y{r) and (j){t,r). We thus use the chain rule: 



dS, 



reg 



du! 



-M 



dr 



— + 

duo duj dy 



dy d£ dy' dC 



+ 



+ 



dC 



du dy' du 



+ 



d4>' dC 
duj dS' 



(20) 



-27- 



Notice in particular that in the first term the acts only on the explicit uj dependence in 

y and 



C We then use the fact that mixed partial derivatives commute to write 
similarly for </>, and then integrate by parts to find 



d_d_ 

dr du> 



reg 



-M J dr 

dy dC 
duo dy' 



+ 



duo 
A 

+ 



dy 



+ 



dC 

dy 



d_dC 

dr dy' 



dy 



duj 



+ 



dC 



d_ dC^ 

dr d(j)' 



duj 



dC 



duj < 



(21) 



Clearly the coefficients of the and tj£ terms inside the integral vanish on-shell. 



Turning to the boundary terms, we start with those for y(r). Notice first that the 



contribution from the r = 
the square root in eq. M 



endpoint vanishes because ^ 



due to the factor outside 



The contribution from the r = A endpoint also vanishes. To see 



this we must use the fact that in the SYM theory m and /j, are independent parameters, so 
that in the supergravity theory — 



can then see that 



dy 

du> 



0. From the asymptotic form of y{r) in eq. (3.10) we 



at r 



A is order That, combined with the fact that in eq. 



(11) is order one (in the A counting), indicates that the contribution from the r = A endpoint 



oo. 



is order ^^^^ and hence vanishes as A 

For the (j) boundary terms we first identify ^ = c. Recalling that our ansatz is (t){t, r) = 
iot + f{r), the factor ^ will give us a term which is simply t, which vanishes since ct is 
independent of r, so (ct)|Q = 0. As for the r dependence in (p, the contribution from the 
r = A endpoint vanishes, as we can see from the asymptotic form of (p in eq. (|1|): the 
leading term is order and hence vanishes as A ^ oo. All that remains is the contribution 
from the r = endpoint. 

We are only interested in cq = solutions, however, for which c = anyway, so that 
the (p boundary terms do not contribute at all. In those cases, the only contribution to (J*) 
comes from the first term under the integral. 

As for the counterterms, Li and L3 contribute nothing, while from the explicit form of 
L2 in eq. (S) we have that 



— -4^^o-(log(cg) 



log(A2) + 1) + O 



log A 



Notice that this vanishes for the cq = solutions that we want. 

To summarize: for the solutions we want, which have c = and cq 
solution for (/>'(r) in eq. (^I^), (J*) is given by 



0, and using the 



(J*) = -Af 



dr— = -M 
duj 



j dr r^ Vl + : 



,/2 



1 + 



1 



■^,2 y 



(jj2 _j_ ^2^2 



(22) 



Notice also that very similar arguments apply for the magnetization [16], given in the 
SYM theory by — The biggest changes are that u; — > 5, the ^ term in y(r)'s asymptotic 
behavior is relevant to show that the y boundary term vanishes at r = A rather than the ^^^^ 
term, and now counterterm L3 contributes rather than L2. 
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